Quasiparticle transport and localization in high-T c superconductors 



T. Senthil, Matthew P.A. Fisher, Leon BalentsJH Chetan Nayaki 
Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106-1^030 

(February 1, 2008) 

We present a theory of the effects of impurity scattering in d x 2_ y 2 superconductors and their 
quantum disordered counterparts, based on a non-linear sigma model formulation. We show the 
existence, in a quasi-two-dimensional system, of a novel spin-metal phase with a non-zero spin 
diffusion constant at zero temperature. With decreasing inter-layer coupling, the system undergoes 
a quantum phase transition (in a new universality class) to a localized spin-insulator. Experimental 
implications for spin and thermal transport in the high-temperature superconductors are discussed. 



Over the last few years, experiments || have convinc- 
ingly established that the superconducting state of the 
hole-doped cuprate materials is characterized by spin 
singlet d x 2_ y 2 pairing. In such a superconductor, the 
gap vanishes at four points on the (two-dimensional) 
Fermi surface. The quasiparticle excitations at these 
"nodes" have a linear dispersion, and an associated den- 
sity of states that vanishes linearly on approaching the 
Fermi surface. This leads to power law dependences of 
various physical quantities on temperature. Impurity 
scattering is expected to strongly modify these proper- 
ties. Experimentally, the power law temperature depen- 
dences are rounded off, apparently approaching constant, 
temperature-independent behavior at the lowest temper- 
atures. This fact is well reproduced by approximate, 
self-consistent treatments of impurity scattering which 
show that a constant finite density of states is generated 
at the Fermi energy for any arbitrarily weak impurity 
strength QB]. Quasiparticle transport properties have 
also been investigated d]?]] theoretically with such self- 
consistent approximations with some phenomenological 
success. Going further, Lee Q has suggested, on the ba- 
sis of calculations of the zero frequency microwave con- 
ductivity, that the quasiparticle eigenstates are strongly 
localized. 

In this paper, we reconsider the effects of disorder on 
the low energy quasiparticles in the d x 2_ y 2 supercon- 
ductor. We point out that the problem of quasiparticle 
transport and localization in a superconductor is con- 
ceptually very different from the more familiar situation 
of non-interacting electrons in a random potential. This 
is because, unlike in a normal metal, the charge of the 
quasiparticles in the superconductor is not a conserved 
quantity. This immediately implies that the quasiparti- 
cle charge in the superconductor cannot be transported 
through diffusion. Indeed the quasiparticle charge den- 
sity is not a hydrodynamic mode in the superconductor. 
However, in a singlet superconductor (and in particular 
in the high-T c superconductors) , the condensate does not 
carry any spin, and consequently the spin of the quasi- 
particles is a good quantum number and is conserved. 
The quasiparticle energy is also conserved. Thus, there 
is the possibility of having spin and energy diffusion with- 
out charge diffusion. These differences in symmetry lead 



to interesting differences between the localization prop- 
erties of quasiparticles in the superconductor, and in the 
normal metal. Such differences have been pointed out 
before Q in the context of the random matrix theory of 
mesoscopic normal/superconducting systems. 

We address quasiparticle transport using a replica field 
theoretic formulation. As expected, the field theory is 
different from that describing Anderson localization in a 
normal metal. The properties of the theory are deter- 
mined by a single coupling constant, which is the dimen- 
sionless spin conductance. This is the physically correct 
quantity whose behavior as a function of system size en- 
ables construction of a scaling theory of localization. By 
analyzing the properties of the field theory, we show the 
existence of a logarithmic "weak localization" correction 
in two dimensions suggesting localization at the largest 
length scales. This correction persists, in part, in the 
presence of an orbital magnetic field (unlike usual An- 
derson localization) or a Zeeman field, but is suppressed 
when both are present. In all cases, however, the quasi- 
particles are generically ultimately localized in two di- 
mensions. Upon inclusion of interlayer coupling, there is 
the interesting possibility of a quantum phase transition 
between an extended spin metal and a localized spin in- 
sulator. The spin metal has diffusive spin correlations, 
a finite spin susceptibility, and an associated finite spin 
conductivity all at zero temperature. We are not aware 
of the existence of such a spin phase in any insulating 
Heisenberg spin model with or without randomness. 

The spin insulator is expected to exhibit local moments 
and spin-glass or random-singlet behavior at very low 
temperatures. The transition between these two phases 
is described by the critical point of the replica field the- 
ory (neglecting quasiparticle interactions), and is a new 
universality class for localization. 

Most of these results also go over unmodified to the 
quantum disordered version of the d x 2_ y 2 superconduc- 
tor - the "nodal liquid" phase that has been analyzed 
recently || as a possible low temperature theory of the 
pseudo-gap regime in the cuprate materials. 

We begin our analysis with the lattice quasiparticle 
Hamiltonian for a singlet superconductor, 
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where i,j are site labels. Using Hermiticity combined 
with spin- rotational and time-reversal invariances, t and 
A may be taken to be real-symmetric matrices. Note 
that the total number of particles is not conserved by this 
Hamiltonian while their total spin is. In the presence of 
impurities, we define U 3 = t° + and A y - = A°- + A} } , 
where t° and A are the Fourier transforms of the ki- 
netic energy (measured from the Fermi energy) and 
the gap function A&, respectively. In the d x i_ y i su- 
perconductor of primary interest, we may take A^ = 
Ao(cos k x — cos k y ). We mention in passing that for weak 
impurities, a continuum limit may be taken, focusing on 
wavevectors near the d-wave nodes. The resulting "dirty 
Dirac" Hamiltonian is similar to various models in the lit- 
erature lie] , but differs from previously studied variants 
in that it contains several random anomalous couplings. 
The Hamiltonian above can be regarded as a lattice reg- 
ularization of this continuum effective field theory for the 
d x 2_ y 2 superconductor, hence our results are quite gen- 
eral and not restricted to a BCS approximation. 

The effect of weak randomness can be analyzed by per- 
turbative renormalization group calculations |H| which 
show that the randomness is a (marginally) relevant per- 
turbation. To make progress then, we employ a field- 
theoretic reformulation of the self-consistent treatment 
adopted in earlier works on dirty d-wave superconduc- 
tivity This begins with the standard coherent-state 
functional integral formulation, in which the electron op- 
erators c, ct are replaced by Grassman fields c, c averaged 
with respect to a statistical weight e~ s , where S is an ac- 
tion. As the randomness is independent of time and Ti. 
is quadratic, different pairs of frequencies (lu, ~oj) decou- 
ple, and it is sufficient for our purposes to focus simply 
on uj = 0. 

Several notational conventions are convenient. We de- 
fine four-component fields ipiaa, with ipn a = Ci a /V2 and 
-0i2 Q = i^ v a nCifj I v2. From this point on we adopt a no- 
tation in which f and a matrices act in the particle/hole 
(a) and spin (a) spaces, respectively. A conjugate field 
is then defined by ip i — (Ci/ , i) T , where C = <j v t v . The 
action in these variables appears non- anomalous, 
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At this stage we have also included an infinitesimal imag- 
inary Zeeman field 77, which acts to generate physical cor- 
relation functions. 

To compute disorder-averaged quantities, we replicate 
the fields 7/7 — ^ ^ p , with /j, = 1 . . . n, so that for n — > 
the statistical weight is normalized for each realization of 
the randomness. Physical quantities can now be simply 
expressed. In particular, the spin susceptibility is xo = 



-(l/^Im^^Vf) (no sums) . Angular brackets denote 
both field-theoretic (ip) and disorder averages. The spin 
diffusion constant, D s , can also be determined from the 
Pij , whose Fourier transform is 



"diffusion propagator" 
P(<l) = Pij exp[« • (a 



Xj)] = 8irxo/(D s q 2 ). One has 
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and no replica sum should be taken. 

The ensemble average over t , A can now be imme- 
diately performed, generating a translationally-invariant 
action with non-trivial quartic couplings between differ- 
ent replicas. A more general analysis |1 1| demonstrates 
that the essential features are captured by uncorrelated 

and Ajj = 

ns. U§ij . 



zero-mean local Gaussian fields t\- = t\5 XJ 



A*i5ij with covariances [t}ti] 
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With this choice, the algebra is particularly simple, and 
the quartic interactions can in turn be decoupled via two 
2n x In Hermitian Hubbard-Stratonovich fields, Q and 
P, acting in the spin and replica spaces (diagonal in the 
particle/hole space). The effective action becomes 

s = E ^ Tr [ww) 2 + ( p (*)) 2 ] + 2 E^ [ W) 

i ij 

-p(i) T y + z W z )s. l3 + t%r z + a°-t*] ^, (3) 

where we have suppressed spin and replica indices. 

A saddle-point (in Q and P) analysis of Eq. || recovers 
the conventional self-consistent approximation. In par- 
ticular, one finds Q — 2ttxo< jZ and P = 0. The constant 
Xo appears as an imaginary self-energy, is (the saddle- 
point approximation to) the physical spin-susceptibility, 
and represents a generation of a non-zero "quasiparticle 
density of states (DOS)" due to disorder. The conclusion 
that xo 7^ is amply supported by experiment, leading 
us to believe that this saddle-point is a physically cor- 
rect starting point. The imaginary self-energy also has a 
complementary interpretation as a finite (inverse) elastic 
scattering time I/r e . For times longer than r e , quasipar- 
ticles no longer move ballistically, and we expect diffusion 
of the conserved energy and spin densities. 

Fluctuations around this saddle-point represent both 
diffusion and corrections to it. Near two spatial dimen- 
sions these fluctuations are captured by a Non-Linear 
Sigma-Model (NLctM) treatment. The crucial ingredi- 
ents are the physical (non-statistical) symmetry prop- 
erties of the Hamiltonian, which determine symmetries 
of the replicated action, Eqs. |^. For the SU(2) and 
time-reversal invariant form chosen, the crucial symme- 
try group is Sp{2n) x Sp(2n). In particular, consider the 
transformation tpi — > Uipi, with U = ^[Ua{^ + t v ) + 
Ub{1 ~ tV )]j with Ua,b 2h x 2n unitary matrices in the 
spin and replica spaces satisfying U\ B cr v U a,b = & v ■ Un- 
der this transformation, the other fields rotate according 
to Q + iP -> U\{Q + iP)U B - For r\ = 0, all such rota- 
tions leave S invariant, while this is true for non-zero r\ 
only when Ub — & z Uac z , hence n breaks the symmetry 
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infinitesimally from Sp(2n) x Sp(2n) to Sp(2n). For a 
single replica, note that as Sp(2) ~ SU(2), one of these 
Sp(2) symmetries is just spin rotation invariance. The 
other Sp(2) symmetry is actually a consequence of time 
reversal invariance, and can be traced to the reality of 
the Hamiltonian H. 

The NLtrM is constructed by considering fluctuating 
Sp(2n) x Sp(2n) rotations of saddle-point solutions that 
are slowly-varying in space. In general, these can be 
shown to take the form of an Sp(2n) matrix U(x), with 
Q(x) + iP{x) = ^\ a z U{x). The form of the action for 
U is determined entirely on symmetry grounds, and is 
verified by a direct calculation |fll|| expanding Q and P 
and integrating out non-critical massive modes. We find 
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where U(x) € Sp(2n). This field theory is known as the 
"principal chiral Sp(2n) model" in the field theory liter- 
ature. In contrast to the conventional sigma models used 
to describe the localization of non-interacting electrons, 
here the field variables live on a group manifold instead of 
a coset space. The Sp(2n) x Sp(2n) symmetry acts on U 
via global left and right multiplication with independent 
Sp(2n) matrices. 

The replica-diagonal self-consistent approximation 
used in other work corresponds to keeping only the con- 
figuration U{x) = 1 in the action. Small quadratic fluc- 
tuations around this solution correspond to diffusion, and 
a direct calculation of P^j in this approximation re- 
lates the coupling-constant to the spin-conductance a s , 
to wit | = \<J S - The derivation of the sigma model 
provides an estimate for the bare coupling constant: 
— = -L v f+ v a w ith vw, the Fermi velocity, and t> A , 
the slope of the d x 2_ y 2 gap linearized near the nodes. 
Note that this is independent of the disorder strength. 
A similar result for the zero frequency microwave con- 
ductance was obtained earlier by Lee Q), in particular 
cr(u! = + ) = ■i 2 -(vF/v & )e 2 /h. The difference in the 
velocity-dependence of the prefactors is conceptually sig- 
nificant: the spin-conductance obeys an Einstein relation 
while the microwave conductance cannot. This distinc- 
tion arises because the quasiparticle charge is not a good 
quantum number. 

Consider separately the orbital and Zeeman cou- 
plings to an applied magnetic field. The orbital field 
breaks time-reversal symmetry but not SU(2), and sim- 
ilar manipulations to those above lead ultimately to a 
Sp(2n)/U{n) NLuM, also distinct from the three conven- 
tional universality classes of dirty metals. The Zeeman 
coupling, by contrast, breaks SU{2) invariance, leaving 
only a U(l) spin-rotation symmetry around the field axis 
(say z). Using the particle/hole transformation, c, — ► ct, 
this U(l) symmetry is easily shown to play the same role 
as does the charge-conservation U(l) in conventional lo- 
calization. Consequently, the Zeeman field leads to the 



usual orthogonal sigma model, and Zeeman and orbital 
effects together drive the system to the unitary univer- 
sality class. 

All these field theories exhibit diffusion on length scales 
of order the elastic mean free path £ e . Beyond this 
scale, quantum interference corrections can play an im- 
portant role. They are determined by the renormaliza- 
tion group (RG) equation for g, which for the Sp(2n) and 
Sp(2n)/U(n) sigma models can be found in Refs. Q: 
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where e = d — 2 and we have set n = 0. The num- 
ber a = 1, ^for the Sp(2n), Sp(2n)/U(n) models respec- 
tively. Eq. H describes the evolution of the physical cou- 
pling (and hence the spin conductance) with length scale 
L, which could be either the system size or an inelastic 
thermal cut-off length at finite temperature. Note that 
in two-dimensions (e = 0), g grows logarithmically with 
L, giving an additive logarithmic reduction of the con- 
ductance and signaling a cross-over to localized behavior 
at long distances. Notice that to this order ("one- loop") 
the leading logarithmic correction is not completely sup- 
pressed by an orbital field. This result is in sharp con- 
trast to conventional weak-localization, but is in agree- 
ment with similar observations made in the context of 
the random matrix theory of systems with these sym- 
metries JlJ,^]. Complete suppression of the logarithmic 
correction occurs only with the introduction of Zeeman 
coupling and subsequent crossover to the unitary NLctM. 
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FIG. 1. Schematic phase diagram of the layered dirty 
d x 2_ y 2 superconductor. 

These perturbative results strongly suggest that in two 
spatial dimensions with weak magnetic fields, the quasi- 
particles are ultimately always localized. A crude esti- 
mate for the localization length may be obtained from the 
one loop perturbation theory to be £ ~ l e e g o where £ e 
is the mean free path and go the bare coupling constant. 
Using the estimate %f ~ 7 in high-T c , we get £ « 10004. 

We note in passing that since II2 ( S y^y^ ) — Z, a. non- 
trivial topological term is allowed for non-zero orbital 
coupling; this suggests the possibility of isolated ex- 
tended quasiparticle states for strong magnetic fields [O . 
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Inclusion of coupling between two-dimensional layers 
(with spacing d) drives the system three dimensional, 
making possible an extended phase where the spins dif- 
fuse at the longest length scales. Based on the quasi- 
2d NLcrM Q , the boundary between 3d spin-metal and 
spin-localized phases occurs when the bare z-axis spin 
conductivity a 7 f ~ d£~ 2 ! e~ 871 : < /g (see Fig. |l]). Given this 
steep curvature of the phase boundary near the origin, 
even a modest interlayer coupling can drive the system 
into the spin metal phase. In zero field, or neglecting 
Zeeman coupling, the transition between the spin metal 
and the spin insulator is in a new universality class. 

Quasiparticle interactions, which we have ignored so 
far, can be shown |Q to lead to the usual Altshuler- 
Aronov singularities for the tunneling density of states for 
the diffusive spin metal. Interaction effects are expected 
to be more crucial in the spin insulator, and ultimately 
should produce a low density (considerably less than, e.g. 
the hole doping) of local magnetic moments which may 
then at low temperature freeze into a spin glass or stay 
paramagnetic in a random-singlet phase with a diverging 
spin susceptibility. 

An important application of the theory outlined here 
is to the quantum disordered d x 2_ y 2 superconductor - 
a novel zero temperature phase that has been proposed 
Q very recently to exist between the antiferromagnetic 
and superconducting regions of the high-T c phase dia- 
gram. The low temperature spin and thermal transport 
properties of the nodal liquid are identical to that of the 
superconductor, and with (weak) disorder, all the results 
mentioned above obtain. One quantitative point is worth 
mentioning: as one moves from the superconductor to- 
wards the antiferromagnet through the nodal liquid, the 
ratio Vf/va decreases monotonically, thereby decreas- 
ing the bare spin conductance. Thus localization effects 
are expected to become more important on going to the 
nodal liquid region. It is interesting that it is precisely 
in this region that experiments find a spin glass phase at 
low temperature. 

We conclude with a few brief suggestions for experi- 
ments, leaving a detailed discussion to Ref. M], Spin 
transport can be probed by NMR techniques. It should 
also be possible to observe the localization physics in 
thermal transport. Ignoring the weak interaction effects, 
we predict that the thermal conductivity k is related to 
the spin conductivity by the Wiedemann-Franz law: 



K/Ta s = 4?r 2 /3. 
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Physically, this follows from the equality of spin and ther- 
mal diffusion constants, the Einstein relation, and the re- 
lation between specific heat and density of states. (The 
Lorenz number differs by a factor of four from the usual 
one as the charge e in the usual formula is replaced by 
spin i in our case.) Within the self-consistent theory this 
gives k — 2k 2 B T(v 2 F + v\) /3vpV£,Ti. In contrast, the mi- 
crowave conductivity does not satisfy an Einstein relation 



and is, in general, not related to the thermal conductivity 
by the Weidemann-Franz law (as can be explicitly seen 
in the self-consistent theory, unless in the limit vp 3> v A 
§)• 

Finally, we expect that localization effects should be 
most pronounced when the nodal anisotropy vf/v & is 
minimized, as is expected to occur on reducing the hole 
concentration within the nodal liquid phase. If signa- 
tures of localization can be observed, it may be useful to 
perturb the system with a Zeeman (i.e. in-plane) field. A 
large enough Zeeman coupling is theoretically expected 
to open the d-wave nodes into Fermi pockets [Q, dra- 
matically increasing the density of states and the bare 
conductance, hence potentially probing some of the lo- 
calization transitions discussed here. 
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